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Abstract 

This paper provides a probabilistic and statistical comparison of the log-GARCH and 
EGARCH models, which both rely on multiplicative volatility dynamics without positivity con- 
straints. We compare the main probabilistic properties (strict stationarity, existence of moments, 
tails) of the EGARCH model, which are already known, with those of an asymmetric version 
of the log-GARCH. The quasi-maximum likelihood estimation of the log-GARCH parameters 
is shown to be strongly consistent and asymptotically normal. Similar estimation results are 
only available for particular EGARCH models, and under much stronger assumptions. The 
comparison is pursued via simulation experiments and estimation on real data. 
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1 Preliminaries 



Since their introduction by Engle (1982) and Bollerslev (1986), GARCH models have attracted much 
attention and have been widely investigated in the literature. Many extensions have been suggested 
and, among them, the EGARCH (Exponential GARCH) introduced and studied by Nelson (1991) 
is very popular. In this model, the log-volatility is expressed as a linear combination of its past 
values and past values of the positive and negative parts of the innovations. Two main reasons 
for the success of this formulation are that (i) it allows for asymmetries in volatility (the so-called 
leverage effect: negative shocks tend to have more impact on volatility than positive shocks of the 
same magnitude), and (ii) it does not impose any positivity restrictions on the volatility coefficients. 

Another class of GARCH-type models, which received less attention, seems to share the same 
characteristics. The log-GARCH(p,q) model has been introduced, in slightly different forms, by 
Geweke (1986), Pantula (1986) and Milh0j (1987). For more recent works on this class of models, 
the reader is referred to Sucarrat and Escribano (2010) and the references therein. The (asymmetric) 
log-GARCH(p, q) model takes the form 

et = crtVt, 

< loga? = w + ELi («i+ 1 { et _ ! >0} +a i _l {et _. <0} )loge?_ i (1.1) 

where at > and (rjt) is a sequence of independent and identically distributed (iid) variables such 
that ErjQ = and ErjQ = 1. The usual symmetric log-GARCH corresponds to the case a + = cx_, 
with ol + = (ai + , . . . , a q+ ) and a_ = (ai_, . . . , a g _). 

Interesting features of the log-GARCH specification are the following. 

(a) Absence of positivity constraints. An advantage of modeling the log-volatility rather 
than the volatility is that the vector 6 = (uj, a + , j3) with j3 = (/3i, . . . , f3 p ) is not a priori subject 
to positivity constraints 1 . This property seems particularly appealing when exogenous variables are 
included in the volatility specification (see Sucarrat and Escribano, 2012). 

(b) Asymmetries. Except when aj+ = Oj_ for all i, positive and negative past values of 

et have different impact on the current log-volatility, hence on the current volatility. However, 

given that loge^ can be positive or negative, the usual leverage effect does not have a simple 

1 However, some desirable properties may determine the sign of coefficients. For instance, the present volatility is 
generally thought of as an increasing function of its past values, which entails f3j > 0. The difference with standard 
GARCH models is that such constraints are not required for the existence of the process and, thus, do not complicate 
estimation procedures. 
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characterization, like Oj + < oij_ say. Other asymmetries could be introduced, for instance by 



is not the case of Model (1.1) except in the symmetric case. 

(c) The volatility is not bounded below. Contrary to standard GARCH models and most 
of their extensions, there is no minimum value for the volatility. The existence of such a bound 
can be problematic because, for instance in a GARCH(1,1), the minimum value is determined 
by the intercept ui. On the other hand, the unconditional variance is proportional to co. Log- 
volatility models allow to disentangle these two properties (minimum value and expected value of 
the volatility). 

(d) Small values can have persistent effects on volatility. In usual GARCH models, 
a large value (in modulus) of the volatility will be followed by other large values (through the 
coefficient (3 in the GARCH(1,1), with standard notation). A sudden rise of returns (in module) 
will also be followed by large volatility values if the coefficient a is not too small. We thus have 
persistence of large returns and volatility. But small returns (in module) and small volatilities are 
not persistent. In a period of large volatility, a sudden drop of the return due to a small innovation, 
will not much alter the subsequent volatilities (because (3 is close to 1 in general). By contrast, as 
will illustrated in the sequel, the log-GARCH provides persistence of large and small values. 

(e) Power-invariance of the volatility specification. An interesting potential property 
of time series models is their stability with respect to certain transformations of the observations. 
Contemporaneous aggregation and temporal aggregation of GARCH models have, in particular, 
been studied by several authors (see Drost and Nijman (1993)). On the other hand, the choice of 
a power-transformation is an issue for the volatility specification. For instance, the volatility can 
be expressed in terms of past squared values (as in the usual GARCH) or in terms of past absolute 
values (as in the symmetric TGARCH) but such specifications are incompatible. On the contrary, 
any power transformation \at\ s (for s ^ 0) of a log-GARCH volatility has a log-GARCH form (with 
the same coefficients in 6, except the intercept u> which is multiplied by s/2). 

The log-GARCH model has apparent similarities with the EGARCH(p, £) model defined by 



under the same assumptions on the sequence (%) as in Model (1.1). These models have in common 
the above properties (a), (b), (c) and (e). Concerning the property in (d), and more generally 
the impact of shocks on the volatility dynamics, Figure 1 illustrates the differences between the two 



replacing oj by Ef=i a i+l{et-i>0} + a i-l{e t _ 4 <o}- The model would thus be stable by scaling, which 




(1.2) 
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models (and also with the standard GARCH). The model coefficients have been chosen to ensure 
the same long-term variances when the squared innovations are equal to 1. Interestingly, a shock 
close to zero has a very persistent impact on the log-GARCH volatility, contrary to the EGARCH 
and GARCH volatilities. However, a sequence of shocks close to zero does impact the EGARCH 
volatility (but not the GARCH volatility). 










-e- GARCH 
-B- EGARCH 

Log-GARCH 
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Figure 1: Curves of the impact of shocks on volatility. The top- left graph shows that a large shock 
has a (relatively) small impact on the log-GARCH, a large but transitory effect on the EGARCH, 
and a large and very persistent effect on the classical GARCH volatility. The top-right graph shows 
the effect of a sequence of tiny innovations: for the log-GARCH, contrary to the GARCH and 
EGARCH, the effect is persistent. The bottom graph shows that even one tiny innovation causes 
this persistence of small volatilities for the log-GARCH. 

This article provides a probability and statistical study of the log-GARCH, together with a 
comparison with the EGARCH. While the stationarity properties of the EGARCH are well-known, 
those of the asymmetric log- G ARCH (p, g) model (1.1) have not yet been established, to our knowl- 
edge. As for the quasi-maximum likelihood estimator (QMLE), the consistency and asymptotic 
normality have only been proved in particular cases and under cumbersome assumptions for the 
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EGARCH, and have not yet been established for the log-GARCH. Finally, it seems important to 
compare the two classes of models on typical financial series. The distinctive features of the two 
models may render one or the other formulation more adequate for certain types of series. 

The remainder of the paper is organized as follows. Section 2 studies the existence of a solution 
to Model (1.1). Conditions for the existence of log- moments are derived, and we characterize the 
leverage effect. Section 3 is devoted to the tail properties of the solution. In Section 4, the strong 
consistency and the asymptotic normality of the QMLE are established under mild conditions. 
Section 6 presents some numerical applications on simulated and real data. Proofs are collected in 
Section 7. Section 8 concludes. 

2 Stationarity, moments and asymmetries of the log-GARCH 

We start by studying the existence of solutions to Model (1.1). 
2.1 Strict stationarity 

Let 0^ denote a /c-dimensional vector of zeroes, and let If. denote the /c-dimensional identity matrix. 
Introducing the vectors 

e t,g = (1{«>o} logef, • • • , l{ et _ 9+1>0 } loge 2 t _ q+1 )' G R q , 

e t, q = (l{ et <o}loge t 2 ,...,l {€( _ g+1< o}loge t 2 _ g+1 )'eR 9 , 

*t = (4, q , e t, q , log of, ... , log al p+1 )' G R 2 *+* 

b t = ((w + log?? 2 )l {j , t>0} ,0;_ 1 ,(6j + log 77 2 )l w<0} ,0 , 9 _ 1 ,w,0;„ 1 ) / G 



and the matrix 



C t 



\ 

we rewrite Model (1.1) in matrix form as 



^{rit>0} a + 


1 {??t>0}«- 




Iq-1 Og-1 


°{q-l)xq 


°(?-l)x P 


l{r]t<0} a + 


1 {r lt <0} a ~ 


1 {Vt<0}P 


°(q-l)xq 


Iq-1 Og-1 


°(q-l)xp 


a + 


a 


P 


°(P-I)xq 


°(p-l)xq 





\ 



(2.1) 



z t = C t z t -i + b t . 



(2.2) 
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We have implicitly assumed p > 1 and q > 1 to write Ct and b t , but obvious changes of notation 
can be employed when p < 1 or q < 1. Let 7(C) be the top Lyapunov exponent of the sequence 
C = {Ct, t G Z}, 

7 (C) = lim -E (log ||CtC t _! . . -dH) = inf - £(log \\C t C^i . . . d||). 

t— »oo t t>l t 

The choice of the norm is obviously unimportant for the value of the top Lyapunov exponent. 
However, in the sequel, the matrix norm will be assumed to be multiplicative. Bougerol and Picard 
(1992a) showed that if an equation of the form (2.2) with iid coefficients (C t ,b t ) is irreducible 2 
and if £'log + ||Co|| and Elog + ||&o|| are finite, 7(C) < is the necessary and sufficient condition 
for the existence of a stationary solution to (2.2). Bougerol and Picard (1992b) showed that, for 
the univariate GARCH(p, q) model, there exists a representation of the form (2.2) with positive 
coefficients, and for which the necessary and sufficient condition for the existence of a stationary 
GARCH model is 7(C) < 0. The result can be extended to more general classes of GARCH models 
(see e.g. Francq and Zako'ian, 2010a). The problem is more delicate with the log-GARCH because 
the coefficients of (2.2) are not constrained to be positive. The following result and Remark 2.1 
below show that 7(C) < is only sufficient. The condition is however necessary under the mild 
additional assumption that (2.2) is irreducible. 

Theorem 2.1. Assume that E log + | log rj^ \ < 00. A sufficient condition for the existence of a 
strictly stationary solution to the log-GARCH model (1.1) is 7(C) < 0. When 7(C) < there exists 
only one stationary solution, which is non anticipative and ergodic. 

Example 2.1 (The log-GARCH(l,l) case). In the case p = q = 1, omitting subscripts, we 
have 

P ) II ( a + 1 {*>0} + «-1{i»<0} + 0) ■ 



CtCt~\ ■ ■ ■ C\ 



{Vt>0} 
1 {'?t<0} 

1 



t-i 



Assume that £'log + | log 770 1 < 00, which entails P(r]Q = 0) = 0. Thus, 



7(C) = £log|aq-l {??0> o} + a-l{ vo<0} +fi\= log |/3 + a+\ a \/3 + a_ 



il-a 



where a = -P(ryo > 0). The condition |a + + /3| a |a_ + fi\ l a < 1 thus guarantees the existence of a 
stationary solution to the log-GARCH(l,l) model. 



2 See their Definition 2.3. 
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Example 2.2 (The symmetric case). In the case a + = a_ = a, one can see directly from (1.1) 
that log of satisfies an ARMA-type equation of the form 

1 1 - (a; + ft) B l 1 log a 2 t =c + Y, 

I i=l J i=l 

where B denotes the backshift operator, vt = logT/ 2 , r = m&x{p,q}, on = for i > q and ft = for 
i > p. This equation is a standard ARMA(r, g) equation under the moment condition E^ogrjf) 2 < 
oo, but this assumption is not needed. It is well known that this equation admits a non degenerated 
and non anticipative stationary solution if and only if the roots of the AR polynomial lie outside 
the unit circle. 

We now show that this condition is equivalent to the condition 7(C) < in the case q = 1. 
Let P be the permutation matrix obtained by permuting the first and second rows of 12+p- Note 
that Ct = C + l{^ t> o} + C _ l/ r?t< o} with C~ = PC + . Since a + = a_, we have C + P = C + . Thus 
C + C = C + PC + = C + C + and ||C t • • • Ci|| = ||(C + )'||. It follows that 7(C) = logp(C + ). In 
view of the companion form of C + , it can be seen that the condition p(C + ) < 1 is equivalent to 
the condition z — Yll=i( a i + ft) z * = =^ |z| > 1. 

Remark 2.1 (The condition 7(C) < is not necessary). Assume for instance that p = q = 1 
and a+ = a_ = a. In that case 7(C) < is equivalent to \a + j3\ < 1. In addition, assume 
that t/q = 1 a.s. Then, when a + /3 7^ 1, there exists a stationary solution to (1.1) defined by 
et = exp(c/2)rjt, with c = w/(l — a — fi). 

2.2 Existence of log-moments 

It is well known that for GARCH-type models, the strict stationarity condition entails the existence 
of a moment of order s > for | et \ . The following Lemma shows that this is also the case for | log e\ \ 
in the log-GARCH model, when the condition E log + | log t/q | < 00 of Theorem 2.1 is slightly 
reinforced. 

Proposition 2.1 (Existence of a fractional log-moment). Assume that 7(C) < and that 
E\ log t/qI 50 < 00 for some sq > 0. Let et be the strict stationary solution of (1.1). There exists 
s > such that E\ log ef | s < 00 and E\ logoff < 00. 

In order to give conditions for the existence of higher-order moments, we introduce some addi- 
tional notation. Let ej be the i-th column of I r , let o"t r = (log of , . . . , log o^_ r+1 )', r = max{p, q}, 
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and let the companion matrix 

Mi(^-l) 




(2.3) 



where m(rjt) = aj+lj^x)} +aj_l{ r?t< o} + Pi with the convention oti+ = a%- = for i > p and fii = 
for i > q. We have the Markovian representation 

(T t ,r = A t (Tt-l,r + U t , (2.4) 

where = itjei, with 

9 

« t = w + ^ (a i+ l {% _ i>0 } + a>i^l {m _ t<0} ) log7/ 2 „ r 

i=l 



The Kronecker matrix product is denoted by ®, and the spectral radius of a square matrix Af is 
denoted by p{M). For any (random) vector or matrix M, let Abs(AT) be the matrix, of same size 
as M, whose elements are the absolute values of the corresponding elements of M. 

Proposition 2.2 (Existence of log-moments). Letm be a positive integer. Assume thatj(C) < 
and that E\ \ogrft\ m < oo. Let = £Abs(Ai)® m where A t is defined by (2.3). 

• Ifm = l or r = 1, then p(A^ m ^) < 1 implies that the strict stationary solution of (1.1) is such 
that E\ loge 2 | m < oo and E\ log erf \ m < oo. 

• If p(C im) ) < 1, then E\ loge?| m < oo and E\ log erf | m < oo. 

Remark 2.2 (A sufficient condition for the existence of any log-moment). Let A^ 00 -* = 
ess sup Abs(Ai) be the essential supremum of Abs(Ai) term by term. Then, it follows from (7.2) 
in the proof that componentwise we have 

oo 

AbsK r ) < ^(A^Abs^). (2.5) 
e=o 

Therefore, the condition 

p(A(°°)) < 1 (2.6) 

ensures the existence of E\ loge 2 | m at any order m, provided 7(C) < and E\ logr^!" 1 < 00. Now 
in view of the companion form of the matrix A^°°> (see e.g. Corollary 2.2 in Francq and Zako'ian, 
2010a), (2.6) holds if and only if 

r 

^2\a i+ + pi\V\ ai -+ fo\ <1. (2.7) 

8=1 
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Example 2.3 (Log-GARCH(l,l) continued). In the case p = q = 1, we have 
A t = a + l {vt _ 1>0} + a_l {vt _ 1<0} + f3 and = E (| A x \) m . 

The conditions E\ log?7g| m < oo and 



in 



m I (a\a + \ k + {l-a)\a-\ k )\P\ m - k <l 
thus entail E\ \oge^\ m < oo for the log-GARCH(l,l) model. 

Example 2.4 (Symmetric case continued). When a + = a_ = a, the matrix A t is no more 
random: 

AM = A« = Abs(Ai) = 1 1 ' ' 

\ Ir-l Or-1 

In view of the companion form of this matrix we have p ( A^ ) < 1 if and only if 



t=i 

The previous condition ensures E\ loge^| m < oo for all m such that E\ logr/pl™ < oo. 
2.3 Leverage effect 

A well-known stylized fact of financial markets is that negative shocks on the returns impact future 
volatilities more importantly than positive shocks of the same magnitude. In the Log-GARCH(l,l) 
model with a_ > max{0, the usual leverage effect holds for large shocks (at least larger than 1) 
but is reversed for small ones. A measure of the average leverage effect can be defined through the 
covariance between i] t _i and the current log- volatility. We restrict our study to the case p = q = 1, 
omitting subscripts to simplify notation. 

Proposition 2.3 (Leverage effect in the log-GARCH(l,l) model). Consider the log- 
GARCH(1,1) model under (2.6). Assume that the innovations rjt are symmetrically distributed, 
E[\ logr/ol 2 ] < oo and + + |a_|) < 1. Then 

cov{n t ^\oga 2 t ) = -(a + - a_) {E{\i] \)t + E(\ Vo \ log rfc)} , (2.8) 

where 

2 w + H a + -a-Wog^o) 
r = Eloga t = — — . 

1 - p - i(a+ + a_) 
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Thus, if the left hand side of (2.8) is negative the leverage effect is present: past negative innova- 
tions tend to increase the log- volatility, and hence the volatility, more than past positive innovations. 
However, the sign of the covariance is more complicated to determine than for other asymmetric 
models: it depends on all the GARCH coefficients, but also on the properties of the innovations 
distribution. Interestingly, the leverage effect may hold with q + > a.— . Simple calculation shows 
that for the EGARCH(1,1) model, cov(r? t _i, log of) = 71. 

3 Tail properties of the log-GARCH 

In this section, we investigate differences between the EGARCH and the log-GARCH in terms of 
tail properties. 

3.1 Existence of moments 

We start by characterizing the existence of moments for the log-GARCH. The following result is an 
extension of Theorem 1 in Bauwens et al, 2008, to the asymmetric case (see also Theorem 2 in He 
et al, 2002 for the symmetric case with p = q = 1). 

Proposition 3.1 (Existence of moments). Assume that 7(C) < and that p ^A^ 00 - 1 ^ < 1. 
Letting A = maxi<j< 9 {|ai+| V X^>o IK^ 00 ^!! < °°> assume that for some s > 

E exp js^A Vl)|logr?o|j < 00, (3.1) 

then the solution of the log-GARCH(p,q) model satisfies E\eo\ 2s < 00. 

Remark 3.1. In the case p = q = 1, condition (3.1) has a simpler form: 



E 



explsf ; — v I" 1 I v A\ logroll 

P l U-|ai + + /3i| Vai_+ft| /' °'J 



< 00. 



The following result provides a sufficient condition for the Cramer's type condition (3.1). 

Proposition 3.2. If E{\t]q\ s ) < 00 for some s > and i]o admits a density f around such that 
/(y _1 ) = o(|y | 5 ) for 5 < 1 when \y\ — > 00 then E exp(si| log t/q | ) < 00 for some s\ > 0. 

In the case p = q = 1, a very simple moment condition is given by the following result. 

Proposition 3.3 (Moment condition for the log-GARCH(l,l) model). Consider the log- 
GARCH(1,1) model. Assume that log" 1- 1 log 7/q | < 00 and E\t]q\ 2s < 00 for s > 0. Assume 
Pi + cti + € (0, l),/3i + ol\- € (0, 1) and cti + A ai_ > 0. Then ctq and eo have finite moments of 
order s/{a\ + V ai_) and 2s/{a\ + V ai_ V 1) respectively. 
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For an explicit expression of the unconditional moments in the case of symmetric log- 
GARCH(p, q) models, we refer the reader to Bauwens et al. (2008). 

3.2 Regular variation of the log-GARCH(l,l) 

Under the assumptions of Proposition 2.3 we have an explicit expression of the stationary solution. 
Thus it is possible to establish the regular variation properties of the log-GARCH model. Recall 
that L is a slowly varying function iff L(xy)/L(x) — > 1 as x — > oo for any y > 0. A random variable 
X is said to be regularly varying of index s > if there exists a slowly varying function L and 
r € [0, 1] such that 

P(X > x) ~ tx~ s L{x) and P(X < -x) ~ (1 - t)x~ s L(x) x -> +oo. 

The following proposition asserts the regular variation properties of the stationary solution of the 
log-GARCH(l,l) model. 

Proposition 3.4 (Regular variation of the log-GARCH(l,l) model). Consider the log-GARCH(l,l) 
model. Assume that Elog + | log | < oo that rjo is regularly varying with index 2s' > 0. Assume 
Pi + «i+ € (0, l),/?i + OL\- G (0, 1) and ct\ + A a.\- > 0. Then o~q and eo are regularly varying with 
index s'/(«x+ V a i-) an d 2s / /(ai+ V a±- V 1) respectively. 

The square root of the volatility, ctq, thus have heavier tails than the innovations when 
OJ1+ V ai_ > 1. Similarly, in the EGARCH(1,1) model the observations can have a much heavier 
tail than the innovations. Moreover, when the innovations are light tailed distributed (for instance 
exponentially distributed), the EGARCH can exhibit regular variation properties. It is not the 
case for the log-GARCH(l,l) model. 

In this context of heavy tail, a natural way to deal with the dependence structure is to study the 
multivariate regular variation of a trajectory. As the innovations are independent, the dependence 
structure can only come from the volatility process. However, it is also independent in the extremes. 
The following is a straightforward application of Lemma 3.4 of Mikosch and Rezapour (2012). 

Proposition 3.5 (Multivariate regular variation of the log-GARCH(l,l) model). Assume the 
conditions of Proposition 3.4 satisfied. Then the sequence (of) is regularly varying with index 
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s' I (cti + V ct\J). The limit measure of the vector = (af, . . . 
relation on the Borel a -field of (R U {-hooD^/jO^} 



, a^)' is given by the following limiting 



P(a 2 > x) 



S' 




s7(ai+Vai_)-l 



l{ yeie .}dy, x ->■ oo. 



where is i/ie i-t/i umi vector in R and the convergence holds vaguely. 

As for the innovations, the limiting measure above is concentrated on the axes. Thus it is also 
the case for the log-GARCH(l,l) process and its extremes values do not cluster. It is a drawback 
for modeling stock returns when clusters of volatilities are stylized facts. This lack of clustering is 
also observed for the EGARCH(1,1) model in Mikosch and Rezapour (2012), in contrast with the 
GARCH(1,1) model, see Mikosch and Starica (2000). 

We have seen that the extremal behavior of the log-GARCH model is similar to that of the 
EGARCH(1,1) model, except that it has a much lighter tail for the same innovations. It is the 
contrary for the extremely small values. Let us consider the inverse (ej -1 ) of the solution (e^) of 
a log-GARCH model. Then, if the variance of (rj^ 1 ) exists and E{j]^ 1 ) = 0, the process (e^ 1 ) 
satisfies a log-GARCH model. Under the assumptions of Proposition 3.4, the extremal behavior of 
6q 1 is driven by the innovation rj^ 1 which, in the gaussian case, is inverse-gamma distributed with 
parameter (1/2,1/2). Thus eQ 1 is regularly varying with index r = 1/2, implying that extremely 
small values of eo are likely to occur. The time to return to the stationary regime from extremely 
small values is much longer than from extremely large values. 

If now we consider the case of an EGARCH model, the invertibility condition obtained by 
Wintenberger and Cai (2011) imposes that the volatility <Tq be bounded from below. Thus, its 
inverse is bounded from above and the extremal behavior of the return eo is the same as that of 
the innovation r/o- As the innovations have lighter tails than the returns, extremely small values 
are not observed. This explains why the time to return to the stationary regime from extremely 
small values is much shorter than from extremely large values, see Figure 1 (the same reasoning 
also holds for GARCH(1,1) model). 



We now consider the statistical inference. Let ei, . . . , c n be observations of the stationary solution 
of (1-1), where 6 is equal to an unknown value 6q belonging to some parameter space C M. d , with 



4 Estimating the log-GARCH by QML 
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d = 2q + p + 1. A QMLE of 0q is defined as any measurable solution n of 



n = arg mm Q n (6), (4.1) 



with 



Q n (0) = n- 1 V 4(0), * + lo g 3?(fl), 

where ro is a fixed integer and log <?( (0) is recursively defined, for t = 1, 2, . . . , n, by 

9 p 
log of (0) = w + ^ (a i+ l {et _ i>0 } + ai-l {e( _ j<0 }) loge?_i + log 5^(0), 

i=i j=i 

using positive initial values for Cq, . . . , ef , ?q(^)' ■ • • > i ^i-p(^)- 

Remark 4.1 (On the choice of the initial values). The initial values can be arbitrary positive 
numbers, for instance eo = ■ ■ ■ = £\-q = <?o(0) = • • • = <7i_ p (0) = v2 (for daily returns of stock 
market indices, in percentage, the empirical variance is often close to 2). They can also depend on 
the parameter, for instance eo = ■ ■ ■ = €-\- q = &o(Q) = • • • = <7i_ p (0) = exp(w/2). It is also possible 
to take initial values depending on the observations, for instance eo = • • • = ei_ g = <7 O (0) = • • • = 



<7j_p(0) = ^n" 1 J2t=i ef • ^ will be shown in the sequel that the choice of ro and of the initial 
values is unimportant for the asymptotic behavior of the QMLE, provided ro is fixed and there 
exists a real random variable K independent of n such that 

sup | log of (0) - loga t 2 (0)| < K, a.s. for t = q -p+ 1, . . . ,q, (4.2) 
eee 

where a\ (0) is defined by (7.3) below. These conditions are supposed to hold in the sequel. 

Remark 4.2 (Our choice of the initial values). Even if the initial values do not affect the 
asymptotic behavior of the QMLE, the finite sample value of 6 n is however quite sensitive to these 
values. Based on simulation experiments and on illustrations on real data, for estimating log- 
GARCH(1,1) models we used the empirical variance of the first 5 values (i.e. a week for daily data) 
as proxi for the unknown value of a\. These initial values allow to compute Sf for t > 2. In order 
to attenuate the influence of the initial value without loosing too much data, we chose ro = 10. 

Remark 4.3 (The empirical treatment of null returns). Under the assumptions of Theo- 
rem 2.1, almost surely ef 7^ 0. However, it may happen that some observations are equal to zero or 
are so close to zero that n cannot be computed (the computation of the logef's being required). 
To solve this potential problem, we imposed a lower bound for the |et|'s. We took the lower bound 
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10 8 , which is well inferior to a beep point, and we checked that nothing was changed in the nu- 
merical illustrations presented here when this lower bound was multiplied or divided by a factor of 
100. 

We now need to introduce some notation. For any 6 £ 0, let the polynomials Ag (z) = 
Yli=i a i,+ zt ' -A.q(z) = Yli=i a i- zl an d Be( z ) = 1 — Yl P j=ifij z ^ ■ By convention, Ag(z) = and 
Ag (z) = if q = 0, and Bg(z) = 1 if p = 0. We also write C(9q) instead of C to emphasize that 
the unknown parameter is 6q. The following assumptions are used to show the strong consistency 
of the QMLE. 

Al: 9q 6 and is compact. 

A2: 7 {C(6>o)}<0 and V# 6 0, \Bg(z)\ = => |jzf| > 1. 

A3: The support of 770 contains at least two positive values and two negative values, Erfc = 1 
and E\ logr/Q| So < 00 for some sq > 0. 

A4: If p > 0, Ag Q (z) and have no common root with Bq {z). Moreover Ag (l)+^4g o (l) 7^ 

and |a g+| + |a g+| + |/3 0p | / 0. 

A5: E I log | < 00. 

Assumptions Al, A2 and A4 are similar to those required for the consistency of the QMLE in 
standard GARCH models (see Berkes et al. 2003, Francq and Zakoian, 2004). Assumption A3 
precludes a mass at zero for the innovation, and, for identifiability reasons, imposes non degeneracy 
of the positive and negative parts of ijq. Note that, for other GARCH-type models, the absence of 
a lower bound for the volatility can entail inconsistency of the (Q)MLE (see Francq and Zakoian 
(2010b) for a study of a finite-order version of the LARCH (00) model introduced by Robinson 
(1991)). This is not the case for the log-GARCH under A5. Note that this assumption can be 
replaced by the sufficient conditions given in Proposition 2.2 (see also Examples 2.3 and 2.4). 

Theorem 4.1 (Strong consistency of the QMLE). Let (0 n ) be a sequence of QMLE satis- 
fying (4-1), where the et's follow the asymmetric log-GARCH model of parameter Oq. Under the 
assumptions (4-2) and A1-A5, almost surely 6 n — > 6q as n — >■ 00. 

Let us now study the asymptotic normality of the QMLE. We need the classical additional 
assumption: 
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o 

A6: 6 £6 and k a := E{t)q) < oo. 

Because the volatility of is not bounded away from 0, we also need the following non classical 
assumption. 

A7: There exists s\ > such that Eexp(s\\ log ?7q | ) < oo, and p(A^°° s> ) < 1. 

The Cramer condition on | logr/pl in A7 is verified if r\t admits a density / around that does not 
explode too fast (see Proposition 3.2). 

Let VQ = (ViQ, . . . , V d Q)' and MQ = QHi.Q', . . . ,M d ,Q')' be the vector and matrix of the 
first-order and second-order partial derivatives of a function Q : — > R. 

Theorem 4.2 (Asymptotic normality of the QMLE). Under the assumptions of Theo- 
rem 4-1 an d A6-A7, we have ^fn{Q n — 6q) -4- jV(0, («4 — 1)J~ 1 ) as n — > oo, where J = 
E[\7 log o~t(Oo)X7 log o~t(0o)'] is a positive definite matrix and -4 denotes convergence in distribu- 
tion. 

It is worth noting that for the general EGARCH model, no similar results, establishing the 
consistency and the asymptotic normality, exist. See however Wintenberger and Cai (2011) for the 
EGARCH(1,1). The difficulty with the EGARCH is to invert the volatility, that is to write of (0) 
as a well-defined function of the past observables. In the log-GARCH model, invertibility reduces 
to the standard assumption on Be given in A2. 

5 Asymmetric log-ACD model for duration data 

The dynamics of duration between stock price changes has attracted much attention in the econo- 
metrics literature. Engle and Russel (1997) proposed the Autoregressive Conditional Duration 
(ACD) model, which assumes that the duration between price changes has the dynamics of the 
square of a GARCH. Bauwens and Giot (2000 and 2003) introduced logarithmic versions of the 
ACD, that do not constrain the sign of the coefficients (see also Bauwens, Giot, Grammig and 
Veredas (2004) and Allen, Chan, McAleer and Peiris (2008)). The asymmetric ACD of Bauwens 
and Giot (2003) applies to pairs of observation (xi,yi), where x% is the duration between two changes 
of the bid-ask quotes posted by a market maker and yi is a variable indicating the direction of change 
of the mid price defined as the average of the bid and ask prices (yi = 1 if the mid price increased 
over duration Xj, and y i = —\ otherwise). The asymmetric log-ACD proposed by Bauwens and 
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Giot (2003) can be written as 

X% — ipi z% , 

< logV'i = Lo + Y,l =1 {ak+^{y^ k =i} + a k _l {yik= _ 1} )logXi_ k (5.1) 

+ Ej=i/ 3 i 1 °g^-i, 

where (zj) is an iid sequence of positive variables with mean 1 (so that ipi can be interpreted as the 
conditional mean of the duration xi). Note that €t := y/xtUt follows the log-GARCH model (1.1), 
with rjt = y/ztDt- Consequently, the results of the present paper also apply to log-ACD models. In 
particular, the parameters of (5.1) can be estimated by fitting model (1.1) on t t = yfxlyt- 

6 Numerical Applications 

6.1 An application to exchange rates 

We consider returns series of the daily exchange rates of the American Dollar (USD), the Japanese 
Yen (JPY), the British Pound (BGP), the Swiss Franc (CHF) and Canadian Dollar (CAD) with 
respect to the Euro. The observations cover the period from January 5, 1999 to January 18, 2012, 
which corresponds to 3344 observations. The data were obtained from the web site 
http : / /www . ecb . int/stats/ exchange/ eurof xref /html/ index . en . html. 

Table 1 displays the estimated log-GARCH(l,l) and EGARCH(1,1) models for each series. For 
all series, except the CHF, condition (2.7) ensuring the existence of any log-moment is satisfied. 
Most of the estimated models present asymmetries. However, the leverage effect is more visible in 
the EGARCH than in the log-GARCH models. This is particularly apparent for the JPY model 
for which 7 is clearly negative. For all models, the persistence parameter j3 is very high. The last 
column shows that for the USD and the GBP, the log-GARCH has a higher (quasi) log-likelihood 
than the EGARCH. The converse is true for the three other assets. A study of the residuals, not 
reported here, is in accordance with the better fit of one particular model for each series. This study 
confirms that the models do not capture exactly the same empirical properties, and are thus not 
perfectly substitutable. 

6.2 A Monte Carlo experiment 

To evaluate the finite sample performance of the QML for the two models we made the following nu- 
merical experiments. We first simulated the log-GARCH(l,l) model, with n = 3344, r\t ~ A^(0, 1), 
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Table 1: Log-GARCH(l,l) and EGARCH(1,1) models fitted by QMLE on daily returns of exchange 
rates. The estimated standard deviation are displayed into brackets. The larger log-likelihood is 
displayed in bold face. 



Log-GARCH 












Q 


a_|_ 


a_ 


P 


Log-Lik. 


USD 


0.024 (0.005) 


0.027 (0.004) 


0.016 (0.004) 


0.971 (0.005) 


-0.104 


JPY 


0.051 (0.007) 


0.037 (0.006) 


0.042 (0.006) 


0.952 (0.006) 


-0.354 


GBP 


0.032 (0.006) 


0.030 (0.005) 


0.029 (0.005) 


0.964 (0.006) 


0.547 


CHF 


0.057 (0.012) 


0.046 (0.008) 


0.036 (0.007) 


0.954 (0.008) 


1.477 


CAD 


0.021 (0.005) 


0.025 (0.004) 


0.017 (0.004) 


0.969 (0.006) 


-0.170 


EGARCH 












Q 


7 


6 




Log-Lik. 


USD 


-0.202 (0.030) 


-0.015 (0.014) 


0.218 (0.031) 


0.961 (0.010) 


-0.116 


JPY 


-0.152 (0.021) 


-0.061 (0.014) 


0.171 (0.024) 


0.970 (0.006) 


-0.334 


GBP 


-0.447 (0.048) 


-0.029 (0.021) 


0.420 (0.041) 


0.913 (0.017) 


0.503 


CHF 


-0.246 (0.046) 


-0.071 (0.022) 


0.195 (0.035) 


0.962 (0.009) 


1.568 


CAD 


-0.091 (0.017) 


-0.008 (0.010) 


0.103 (0.019) 


0.986 (0.005) 


-0.161 



and a parameter close to those of Table 1, that is 6q = (0.024, 0.027, 0.016, 0.971). Notice that as- 
sumptions A1-A4 required for the consistency are clearly satisfied. Since |/3q| + \ (|«o+| + l a 0-|) < 
1, A5 is also satisfied in view of Example 2.3. The assumptions A6-A7 required for the asymptotic 
normality are also satisfied, noting that |ao+ + A)| V l Q 0- +A)| < 1 an d using Remark 2.2. The first 
part of Table 2 displays the log-G ARCH (1,1) models fitted on these simulations. This table shows 
that the log-GARCH(l,l) is accurately estimated. Note that the estimated models satisfy also the 
assumptions A1-A7 used to show the consistency and asymptotic normality. We also estimated 
EG ARCH (1,1) models on the same simulations. The results are presented in the second part of 
Table 2. As in Wintenberger and Cai (2011), the optimization of the EGARCH(1,1) models has 
been performed under the constraints 5 > I7I and 



E 1 ^ 

i=l 



$\ (7£t-i + «5|e t _i|)exp (-^737j ) 



< 0. 



These constraints guarantee the invertibility of the model, which is essential to obtain a model that 
can be safely used for prediction (see Wintenberger and Cai (2011) for details). Comparing the 
log- likelihood given in the last column of Table 2, one can see that, as expected, the likelihood 
of the log-GARCH model is greater than that of the (misspecified) EGARCH model, for all the 
simulations. 
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Table 2: Log-GARCH(l,l) and EGARCH(1,1) models fitted on 5 simulations of a log-G ARCH (1,1) 
model. The estimated standard deviation are displayed into brackets. The larger log-likelihood is 
displayed in bold face. 



Log- 


GARCH 










Iter 




a_|_ 


s?_ 


P 


Log-Lik. 


1 


0.025 (0.004) 


0.028 (0.004) 


0.018 (0.004) 


0.968 (0.005) 


-0.415 


2 


0.021 (0.003) 


0.023 (0.003) 


0.013 (0.003) 


0.976 (0.004) 


-0.634 


3 


0.026 (0.003) 


0.028 (0.004) 


0.017 (0.003) 


0.969 (0.004) 


-0.754 


4 


0.022 (0.003) 


0.024 (0.004) 


0.018 (0.003) 


0.972 (0.004) 


-0.389 


5 


0.024 (0.003) 


0.028 (0.004) 


0.014 (0.003) 


0.974 (0.003) 


-0.822 


EGARCH 










Iter 


Q 


7 


6 


P 


Log-Lik. 


1 


-0.095 (0.016) 


-0.014 (0.009) 


0.104 (0.017) 


0.976 (0.006) 


-0.424 


2 


-0.127 (0.018) 


0.009 (0.010) 


0.148 (0.021) 


0.976 (0.007) 


-0.645 


3 


-0.147 (0.018) 


0.001 (0.010) 


0.177 (0.022) 


0.971 (0.007) 


-0.770 


4 


-0.136 (0.019) 


-0.012 (0.010) 


0.155 (0.022) 


0.976 (0.007) 


-0.404 


5 


-0.146 (0.019) 


-0.009 (0.010) 


0.177 (0.023) 


0.971 (0.007) 


-0.842 



In a second time, we repeated the same experiments for simulations of an EGARCH(1,1) model 
of parameter (u , 7o, S , /3 ) = (-0.204,-0.012,0.227,0.963). Table 3 is the analog of Table 2 for 
the simulations of this EGARCH model instead of the log-GARCH. The EGARCH are satisfactorily 
estimated, and, once again, the simulated model has a higher likelihood than the misspecified model. 
From this simulation experiment, we draw the conclusion that it makes sense to select the model 
with the higher likelihood, as we did for the series of exchange rates. 

7 Proofs 

7.1 Proof of Theorem 2.1 

Since the random variable ||Co|| is bounded, we have Elog + ||Co|| < oo. The moment condition 
on rjt entails that we also have E\og + ||&o|| < oo. When 7(C) < 0, Cauchy's root test shows that, 
almost surely (a.s.), the series 

00 

z t = b t + (7.1) 

n=0 

converges absolutely for all t and satisfies (2.2). A strictly stationary solution to model (1.1) is then 
obtained as et = exp {5^25+1,*} Vt, where z^t denotes the i-th element of z t . This solution is non 
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Table 3: As Table 2, but for 5 simulations of an EGARCH(1,1) model. 



Log- 


GARCH 










Iter 


Q 




S?_ 


8 


Log-Lik. 


1 


0.039 (0.008) 


0.071 (0.008) 


0.052 (0.007) 


0.874 (0.015) 


-0.350 


2 


0.055 (0.006) 


0.058 (0.007) 


0.052 (0.006) 


0.913 (0.010) 


-0.476 


3 


0.052 (0.008) 


0.070 (0.008) 


0.060 (0.007) 


0.873 (0.015) 


-0.468 


4 


0.051 (0.008) 


0.076 (0.008) 


0.056 (0.007) 


0.878 (0.014) 


-0.416 


5 


0.056 (0.007) 


0.061 (0.007) 


0.060 (0.007) 


0.896 (0.012) 


-0.517 


EGARCH 




8 






Iter 


LO 


7 


P 


Log-Lik. 


1 


-0.220 (0.022) 


-0.024 (0.013) 


0.235 (0.023) 


0.950 (0.010) 


-0.335 


2 


-0.196 (0.020) 


-0.029 (0.012) 


0.219 (0.022) 


0.961 (0.008) 


-0.468 


3 


-0.222 (0.022) 


-0.005 (0.013) 


0.241 (0.024) 


0.947 (0.010) 


-0.448 


4 


-0.227 (0.022) 


-0.025 (0.012) 


0.248 (0.023) 


0.950 (0.010) 


-0.402 


5 


-0.209 (0.021) 


-0.003 (0.012) 


0.234 (0.023) 


0.955 (0.009) 


-0.504 



anticipative and ergodic, as a measurable function of {r] u ,u < t}. 

We now prove that (7.1) is the unique nonanticipative solution of (2.2) when 7(C) < 0. Let 
(zt) be a strictly stationary process satisfying z\ = CtZ^_ 1 + bf. For all N > 0, 

N 

z* = z t (N) + C t ... C t - N z*^ N _ X t z t (N) = b t + ^2 CtCt-x ■ ■ ■ C^ n bt-n-i- 



We then have 



\ z t 



z*t\\ < 



-nPt-n-l 

n=N+l 



n=0 



+ \\c t ...c 



t-N\ 



't-N-l\ 



The first term in the right-hand side tends to a.s. when iV — > 00. The second term tends to in 
probability because 7(C) < entails that \\Ct . . . Cj_tv|| —> a.s. and the distribution of ||^_jy_i|| 
is independent of by stationarity. We have shown that Zt — z\ — > in probability when N — > 00. 



This quantity being independent of N we have z% = z\ a.s. for any t. 



□ 



7.2 Proof of Proposition 2.1 

Let A be a random variable such that X > a.s. and EX r < 00 for some r > 0. If Slog A < 0, 
then there exists s > such that EX S < 1 (see e.g. Lemma 2.3 in Berkes, Horvath and Kokoszka, 
2003). Noting that E \\Ct ■ • • C\\\ < (E ||Ci||)* < 00 for all t, the previous result shows that when 
7(C) < we have E\\Ck • • • C\\\ s < 1 for some s > and some ka > 1. One can always assume 
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that s < 1. In view of (7.1), the c r -inequality and standard arguments (see e.g. Corollary 2.3 in 
Francq and Zako'ian, 2010a) entail that £7||z t || s < oo, provided S||b t || s < oo, which holds true when 
s <■ s o- The conclusion follows. □ 



7.3 Proof of Proposition 2.2 

By (2.4), componentwise we have 

oo I 

Abs(ov) < Abs(u t ) + A t ^Abs(ti t _/_i) J A t/ := JJ Abs(A t _ i ), (7.2) 

where each element of the series is defined a priori in [0, oo]. In view of the form (2.3) of the matrices 
At, each element of 

t 

A t ,eAbs(u t -£~i) = \ut-t-i\ \\ Abs(A t _y)ei 

3=0 

is a sum of products of the form Ilj=o l^ijivt— with < k < £ and < iq < • • • < i\~ < 

£ + 1. To give more detail, consider for instance the case r = 3. We then have 



/ |/tii(?7t_i)||/ii (774-2)1 |t«t-2| + |M2(?7t-2)|K-2| ^ 



A tA Ahs(u t - 



1^1(774-2) 1 1 U t -2 1 

\ K- 2 | / 

Noting that |ut_^_i| is a function of rjt-e-2 and its past values, we obtain EAt t iAbs(ut-2) = 
EAbs(At)EAbs(At-i)EAbs(ut-2)- More generally, it can be shown by induction on I that the 
i-th element of the vector At_i^_iAbs(ttt_^_i) is independent of the i-th element of the first row 
of Abs(Aj). It follows that EA t /Abs(u t -£-i) = EAbs(A t )EA t -i/-iAbs(u t -£-i). The property 
extends to r 7^ 3. Therefore, although the matrices involved in the product A t ^Abs{u t _^_i) are 
not independent (in the case r > 1), we have 

i 

EA t ^Abs{xH-i-\) = JJSAbs(A 4 _ j )^Abs(u t _ £ _i) 

3=0 

A (1) ) EAbs(ux). 



In view of (7.2), the condition p(A^) < 1 then entails that EAbs(at,r) is finite. 

The case r = 1 is treated by noting that A t ^Abs(u t _£_i) is a product of independent random 
variables. 
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To deal with the cases r / 1 and m / 1, we work with (2.2) instead of (2.4). This Markovian 
representation has an higher dimension but involves independent coefficients Cf. Define Cm by 
replacing Af-j by Ct-j in A t £. We then have 

SC® m Abs(6 t _^_ 1 )® m = (c (m) ) £+1 ^Abs(6 1 ) 0m . 

For all m > 1, let ||M|| m = (^UMlp) 1 /" 1 where ||M|| is the sum of the absolute values of 
the elements of the matrix M. Using the elementary relations ||JW||||JV|| = \\M ® iV|| and 
£7||Abs(M)|| = ||£Abs(M)|| for any matrices M and N, the condition p(C (m) ) < 1 entails 
£||C tji ?Abs(b t _^_i)|| m = ||SC^ m Abs(b t _^_i)® m || -> at the exponential rate as £ -> oo, and 
thus 

oo 

||Abs(z t )|| m < ||Abs(6 t )|| m + \\C t/ Ab s (b M )\\ m < oo, 
which allows to conclude. □ 



7.4 Proof of Proposition 2.3 

By the concavity of the logarithm function, the condition \a + + /3||a_ + /3| < 1 is satisfied. By 
Example 2.1 and the symmetry of the distribution of 770, the existence of a strictly stationary 
solution process (e t ) is thus guaranteed. By 2.3, this solution satisfies E\ log£j| < 00. Let 

a t = (a + l {m>0} + a-l ivt<0 y)vt, h = (a + l {vt>0} + a^l {m<0} )r] t log rjf. 

We have 

Ea t = (a+ - a-)E(r] l {vo>0} ), Eb t = (a+ - a-)E(rj log r?ol {j?0>0} ), 
using the symmetry assumption for the second equality. Thus 

cov(77i_i,log(df )) = E[a t - 1 log(at_ 1 ) + b t -i], 
and the conclusion follows. □ 
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7.5 Proof of Proposition 3.1 

By definition, |log(erf)| < l|°"t,r|| = l|Abs(<T ir )||. Then, we have 



Etf\' < E{e X p( S ||AbsK r )||)} = ^ S ll Abs (^)H fc 



k=0 

oo s *"' 



^ fe ||Ab S ( W o)||^{E£olK^ (oo) )l} & 

* E 



k=0 

= £;exp| S ||Abs( Uo )||f;i|(AM^||), 
I e=o J 

where the last inequality comes from (2.5). By definition Uq = (i^O',^)' with 

i 

u = u} + ^2(ai + l v _ i>0 + cxi-V^o) log rfii. 
i=i 

Thus ||Abs(ito)|| < \uq\ < \u\ + maxi<j< g \ai + \ V |aj_| X^f=i I log 1 and it follows that 

£|<7 2 | s < exp| S Mf;||(A^||l{^exp( S A|log % 2 |)} 9 <oo 
under (3.1). □ 

7.6 Proof of Proposition 3.2 

Without loss of generality assume that / exists on [—1,1]. Then there exists M > such that 
/(1/y) < -^lyl* 5 for all y > 1 and we obtain 

^pH^I) < / ex P (2, 1 log(l/ I ))/(^ + /ex P (, 1 log(^) W (x, 

/oo 
y 2( S i-l)+5 dy + £ ;(| r?o |2 S1 )_ 

The upper bound is finite for sufficiently small s± and the result is proved. □ 

7.7 Proof of Proposition 3.3 

We will use Tweedie's (1988) criterion, which we recall for the reader's convenience. Let (Xt) denote 
a temporally homogeneous Markov chain on a state space E, endowed with a cr-field T . 

Lemma 7.1 (adapted from Tweedie (1988), Theorem 1). Suppose [i is a subinvariant measure, that 
is, 



JE 
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and A £ J 7 is such that < f-t(A) < oo. Suppose there exist a nonnegative measurable function V 
on E, and constants K > and c € (0, 1) such that 

i) E[V(X t ) | Xt-i = x] < K, x £ A, 

ii) E[V(X t ) | X t -i = x] < (1 - c)V(x), x e A c . 

Then, jjl is a finite invariant measure for the chain (Xt) and j E Vdfi < oo. 

It should be noted that this criterion does not make any irreducibility assumption. We have Ut = 
u + (ai + l {llt _ 1>0} + ai_l {r)f _ 1<0} )log(?? i 2 _ 1 ), and 

which shows that (cr 2 ) is a temporally homogeneous Markov chain on ]R + *. By Example 2.1, the 
conditions ensuring the existence of a strictly stationary solution are satisfied. The stationary 
distribution thus defines an invariant probability fx for the Markov chain (of). Let A = [0,K] for 
some K > such that [J, (A) > 0. 

The existence of a s-order moment for rfc entails that e Ul admits a moment of order so := 
s/(ai+ V ai_). Let V(x) = 1 + x s ° for x > 0. For any x > and for < c < 1, we have for x £ A c 
with A' sufficiently large, 

E[V^ 2 t) I of-i = x] = 1 + x (/Jl+ai+)s °£[e SoUl l, w> o] + x^ 1+ai - )so E[e SOUl l V0<0 ] 
< (l-c)V(x), 

because /3i + a\+ < 1 and fi\ + ai_ < 1. On the other hand it is clear that E[V(af) \ o~t_\ = x] 
is bounded for x belonging to A. It follows by the above lemma that E^[V(a^)] < oo where the 
expectation is computed with the stationary distribution. □ 

7.8 Proof of Proposition 3.4 

To prove the first assertion, note that if r]Q is regularly varying of index 2s' then r/ 2 is regularly 
varying of index s'. Thus u\ = ui + (cti+l/^ 0> o} + ai_l{,j 0< o}) log(^g) i s such that 

P( e «i > x ) = p( Vo > 0)P {(vo) ai+ > xe~ u | % > 0} 

+P( m < 0)P{( V 2 )^- > xeT" | 770 < 0} . 
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Then e ui is also regularly varying with index s' := s/{a\ + V Note that 

U1 2(/3+ai+) 



P(fff > x) = P(?7o > 0)P (^e u Vo lP+ai+; > x | T/o > 

+P(? ?0 < 0)P (e"Vo (/3+Ql - ) > x | % < o) . 

As e ni admits regular variation of order s' Q , it admits a moment of order Sq(/3 + oi+) < So- Note that 
-E|?7o| 2s < 00 f° r an Y s < s '- An application of Proposition 3.3 thus gives 

for t > small enough. By independence between u\ and ctq conditionally on 770 > 0, we may apply 
a result by Breiman (1965) to conclude that 

P (e^ 2(/3+Ql+) > x I m > 0) ~ E (af +ai+)s '») P(e? > x | % > 0), 

as x —7- 00. Applying the same arguments to P ^e" 1 ^^" 1 " 01- ^ > x | 770 < 0^ we obtain 

P(a 2 > x) ~ P(t?o > 0)£? (ao (/3+ai+)s «) P (e Ul > x | % > 0) 



+P(r ?0 < 0)P ^ 2(/3+Ql - K ) P(e Ul > x I r? < 0) 

and the first assertion follows. The second assertion follows easily by independence of 770 and <7o> 
with respective regularly variation indexes s' and s f . □ 

7.9 Proof of Theorem 4.1 

We will use the following standard result (see e.g. Exercise 2.11 in Francq and Zakoian, 2010a). 

Lemma 7.2. Let (X n ) be a sequence of random variables. If sup n E\X n \ < 00, then almost surely 
n~ 1 X n — > as n — > 00. The almost sure convergence may fail when sup n E\X n \ = 00. If the 
sequence (X n ) is bounded in probability, then n~ 1 X n — > in probability. 

Turning to the proof of Theorem 4.1, first note that A2, A3 and Proposition 2.1 ensure the a.s. 
absolute convergence of the series 

loga 2 (0) := B e \B) jo, + £ {a i+ l {et _ i>0} + a 4 _l {et _ i<0} ) loge 2 ^} • (7.3) 

Let 



i=l 



n 2 

Q n {0) = n- 1 M*)> it{e) = ^— + \oga 2 t {9). 
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Using standard arguments, as in the proof of Theorem 2.1 in Francq and Zakoian (2004) (here- 
after FZ), the consistency is obtained by showing the following intermediate results 

i) lim sup \Q n (6) - Q n {0)\ = a.s.; 

ii) if aj{9) = aj{9 ) a.s. then 9 = 9 ; 

Hi) E£ t (G ) < oo and if 9 ^ 9 , El t {9) > E£ t (0 ); 
iv) any 9 ^ 9q has a neighborhood V(9) such that 
liminf inf Q n (9*) > E£ t (9 ) a.s. 

n->oo 0*£V(0) 

Because of the multiplicative form of the volatility, the step i) is more delicate than in the standard 
GARCH case. In the case p = q = 1, we have 

loga 2 (0) - logSf (0) = {loga?(0) - \ogal(9)} , Vt > 1. 

In the general case, as in FZ, using (4.2) one can show that for almost all trajectories, 



sup log a 2 (0)- log a 2 (0) <Kp\ 
flee 



(7.4) 



where p £ (0, 1) and > 0. When the initial values are chosen as suggested by Remark 4.2, K 
is the realization of a random variable which is measurable with respect to a({r] u ,u < 5}). First 
complete the proof of i) in the case p = q = 1 and q + = a_, for which the notation is more explicit. 
In view of the multiplicative form of the volatility 



i-2 



^2(0) = e ^- 1 logcrf(0)-|-|- e ^{a;+alog6f_ 1 _ i } ) 



(7.5) 



i=0 



we have 



log 



1 *" 2 

-J>{a; + aloge?_i-i} 



8=0 



1 



log e 



-/3*- 1 log <r?(0) 



^-Mogafce) 



Applying Lemma 7.2, the first term of the right-hand side of the equality tends almost surely to 
zero because it is bounded by a variable of the form \Xt\/t, with -E|At | < oo, under A5. The second 
term is equal to 

- t log | {log a\ (9) - log<? 2 (0)} p^e-P- 1 ** , 

where x* is between log<7 2 (0) and logaf(9). This second term thus tends to log |/3| < when 
t — > oo. It follows that 

< Kp\ (7.6) 



sup 
flee 



1 
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where K and p are as in (7.4). Now consider the general case. Iterating (1.1), using the compactness 
of and the second part of A2, we have 

t-l 



loga t 2 (0) = ^q(9)+q + (0)l {£M>o} log ef_ { + q_(0)l{ et _. <O } loge 



2 

t-i 



i=l 



with 



+X> J (0)iog O 3 +1 _ i (0) 

i=i 



supmax{|c i (e)|,|c i+ (e)|,|c i _(fl)|,|c ijl (e)|,...,|c ij p(0)|} < Kp\ pe (0,1). 



We then obtain a multiplicative form for of (0) which generalizes (7.5), and deduce that 



log 



1 



a\ + 02, 



where 



t-l 



a.s. 



a i = -Y^2ci(6) + c i+ (e)l {et i>0} loge 2 ^ + Ci-(0)l {tt i<0} loge 
i=i 

in view of (7.7) and Lemma 7.2, and for x*'s between \oga^ +1 _j(6) and log5 : g +1 _ J (0), 

expi kg ^^.(fl) I -exp J - ^c tJ (0)loga g 2 +w W 



"2 



flog 



3=1 



7 log 



flog 



^ c tJ (0) {log (0) - log a 2 q+1 _ 3 (9)} cxp <j - lo § 



fc=i 



o(l) a.s. 



(7.7) 



using (4.2) and (7.7). Using again (4.2), it follows that lim sup n _ ;>00 ai < logp < 0. We conclude 
that (7.6) holds true in the general case. The proof of i) then follows from (7.4)-(7.6), as in FZ. 
To show it), note that we have 



B e (B) log a\ (0) = u> + A+(B)l {tt>0} log e? + A e (B)l {et<0} log el 



(7. 



If logaf(6) = logcr^(0o) a - s -, by stationarity we have log of (0) = log of (0q) f° r an t, an d thus we 



have almost surely 



U+(B) Aj (B) \ 

< - — - ~ mr" > l{e t >0} io g e 



w 



2 f^(g) ^ (g) ] 

t + \e (i?) b 0o (s)J 



{ ei <o}loge? 



Be (l) B e {iy 
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Denote by R t any random variable which is measurable with respect to a ({i] u ,u < t}). If 

^(B) A + eo (B) A e (B) A ea (B) 

B e (B) T B 6o (B) Bg(B) T B 6o (B)> { ' ; 

there exists a non null (c+, c_)' G M 2 , such that 

c+hvtX*} lo S e t + c -!{r, t <o} !oge 2 + #t-i = a.s. 
This is equivalent to the two equations 

(c+ log 7? 2 + C+ logof + i? t _x) l{r, t>0 } = 

and 

(c_ log rjl + c_ log of + J?t_i) 1{ %<0 } = 0. 

Note that if an equation of the form alogx 2 l{ x> o} + &l{ x >o} = admits two positive solutions 
then a = 0. This result, A3, and the independence between rjt and (af,Rt-i) imply that c+ = 0. 
Similarly we obtain c_ = 0, which leads to a contradiction. We conclude that (7.9) cannot hold 
true, and the conclusion follows from A4. 

Since of (#) is not bounded away from zero, the beginning of the proof of Hi) slightly differs 
from that given by FZ in the standard GARCH case. In view of (7.8), the second part of A2 and 
A5 entail that E\ logof (0)| < oo for all € 9. It follows that El^{6) < oo and E\£ t (G )\ < oo. 

The rest of the proof of Hi), as well as that of iv), are identical to those given in FZ. □ 

7.10 Proof of Theorem 4.2 

A Taylor expansion gives 

V;Q n (0„) - V,Q n (6» ) = M L Q n (e nti )(e n - O ) for all 1 < i < d, 

where the n j's are such that \\0 n ,i ~ 0o|| < ||0n ~~ 00 II- As in Section 5 of Bardet and Wintenberger 
(2009), the asymptotic normality is obtained by showing: 

1. n 1 /2 V Q n (0 o )^AA(O,(K 4 -l)J), 

2. ||EI(3 n (0 n ) — J|| converges a.s. to for any sequence (0 n ) converging a.s. to 6q and J is 
invert ible, 

3. n l / 2 \\ VQ n (O n ) — VQ n (O n )\\ converges a.s. to 0. 
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In order to prove the points 1-3 we will use the following Lemma 

Lemma 7.3. Under the assumptions of Theorem 4-1 and A7, for any m > there exists a neigh- 
borhood V of Oo such that .E[sup v (o 2 /of (0)) m ] < oo and E[sup v | log of (9)\ m ] < oo. 

Proof. We have 

p 

log of (0 O ) - log o 2 (0) = u - co + £ /3,{log of^o) - log 

i=i 

+V eo _ (T t _i )r + -4e o _ (-B)l J7i> o log??? + ^ 0o _ (5)l^ <o logr ?t 2 
with <r t)J . = (log of (0 O ), • • . ,logo 2 _ r+1 (0 o ))', 

= (ai+lj^ 1>0 } + ai_l{^_ 1<0 } a r+ l^ tr>0 y + a r -l{ Vtr<0 y + /3 r ). 

Under A2, we then have 

logo 2 (6> ) - logo 2 (0) = Bg\B) {UJQ-0J + V 9a -B<r t -i,r 

+(^0 o -0( 5 )V>olog% 2 + A go _ (B)l Vt<o log V ?} . 

Under A7 the assumptions of Proposition 3.1 hold. From the proof of that proposition, we thus 
have that £^exp(5||Abs(<T^ r )||) is finite for some 5 > 0. 

Now, note that V#, Aq(1) and Ag(l) are continuous functions of 6. Choosing a sufficiently 
small neighborhood V of #o, one can make supy ||V0 O _0||, sup^ |^4q o _ (1)| and supy |^4q o _ (1)| 
arbitrarily small. Thus £ , [exp(msup v ||V0 () _0<r tjr ||)] and £^[exp(msup v \\(Ag Q _ (B)l rit _ 1> Q + 

■l) ID] are finite for an appropriate choice of V depending on m. We conclude 
that E [exp (msupy |log {o 2 (0o)/°f (0)} |)] < oo an d the first assertion of the lemma is proved. 

Consider now the second assertion. We have 

sup|logo 2 (0)| < |logo f 2 | + sup|log(o 2 (0 o )/o 2 (0))|. 
V V 

We have already shown that the second term admits a finite moment of order m. So does the first 

term, under A7, by Remark 2.2. □ 

Now let us prove that the point 1. follows from the fact that VQ n (0o) is a martingale in L 2 . 

Indeed 

n 

VQ n (0 ) = n- 1 ^ (l-r/t 2 )Vlogo 2 (0 o ). 

t=r +l 
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As rjt is independent of log cr 2 (#o) and Erf = 1 the Central Limit Theorem for martingale differences 

o 

applies whenever Q = (^4 — 1)E(S7 log <7 2 (#o)V log af(0o)') exists. For any SB, the random vector 
Vlogof (0) is the stationary solution of the equation 

p 



Vloga 2 (0) = ^^Vlogaf^C©) + 

3=1 



where <7 2 p (0) = (loga 2 (0), . . . , log a 2 „ p+1 (0))'. 



1 



~t-l,q 



H-l,q 



(7.10) 



Assumption A2 entails that Vlogo" 2 (0) is a linear combination of e + 



t— t—i,q 



and logo- 2 _j(0) 



for i > 1. Lemma 7.3 ensures that, for any m > 0, there exists a neighborhood V of 0q such that 
-Efsupy I log cr 2 _j(0)| m ] < 00. By Remark 2.2, e.f_ iq and e^- admit moments of any order. Thus, 
for any m > there exists V such that ^[supy ||Vlog<T 2 (0)|| m ] < 00. In particular, Vlogo" 2 (#o) 
admits moments of any order. Thus point 1. is proved. 
Turning to point 2., we have 



MQ n (e) 



n 



1 E 1 

t=r +l 



tie), 



where 



m t (e) 



1 



Vlog<7 2 (0)Vloga 2 (0)'. 



(7.1i; 



By Lemma 7.3, the term cj 2 (0o)/°' 2 (^) admits moments of order as large as we need uniformly on a 
well chosen neighborhood V of 0q. Let us prove that it is also the case for H log er 2 (#). Computation 
gives 



v 



Hloga 2 (0) = ^^HlogatjCO) + 

3=1 

From this relation and A2 we obtain 



0(2g+l)xd 

W 2 (0) 



+ 



0(2<j+l)xd 

V'<x 2 (0) 



:ioga 2 (0) 



0(2g+l)xd 

B e {B)^Va*{e) 



+ 



0(2(j+l)xd 

^e(-B) _1 V / cr|_ 1)P (0 



Thus Hlog<7 2 (#) belongs to C(V) and is integrable because we can always choose V such that 
sup v || V / (T 2 _ ljP (0)|| G L m (see the proof of point 1. above). 

An application of the Cauchy-Schwarz inequality in the RHS term of (7.11) yields the integra- 
bility of supy HL£i(0). The first assertion of point 2. is proved by an application of the ergodic 
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theorem on (M£ t (9)) in the Banach space C(V) equipped with the supremum norm: 

sup ||HQ n (0) - E[m (6)] || -»• a.s. 
v 

An application of Theorem 4.1 ensures that 6 n belongs a.s. to V for sufficiently large n. Thus 

\\MQ n (O n ) - E[B£ (e )]\\ < sup \\MQ n (0) - E[B£o(0)}\\ + \\E[m (9 n )} - E\SM (e )]\\ 

v 

converges a.s. to by continuity of — > E\¥LIq{6)] at 6q as a consequence of a dominating 
argument on V. The first assertion of point 2. is proved. 



The matrix J is non invertible iff there exists a non null deterministic vector A = 
(Ao, Ax, ■ ■ ■ , \p-\-2q)' such that A' V log of (Oq) = a.s. If the latter equality holds, (7.10) entails 
(1, ef t q , ej q , O"| )P (0o)) ^ = a.s. By A3 and arguments used for proving the point (ii) in the proof 
of Theorem 4.1, we deduce A p +x = Ap+^+i = 0. For the same reason if Ax = • • • = Aj = then 
Ap+2 = • • • = Ap + x+i = and \ p+q+ 2 = • • • = A p+(?+ x+i = for 1 < i < p A q. Thus as A / 
there exists A& 7^ for 1 < k < q such that A,- = for j < k. Then, log(cr? . ) is a linear combi- 
nation of (log(<7f_ J -))fc<.,-<g, (l{e t _ J ->o}log(ef_ J -))fc<i<p and (l{ et _ j<0 } logie 2 ^))^^. We thus find 
a log-GARCH(j/, q') representation with p' < p and q' < q, in contradiction with A4. Thus J is 
invertible. 

From (7.10) and an equivalent representation for Vlog<7j(0), we have 

p 

Vlog^ 2 (0)- Vlog5f(0) = ^(Vlo gO - t %.(0)-Vlog ^(0)) 




- ?ti lP (») 

where <T 2 p is defined as c 2 p . Thus, there exist continuous functions di and dt.i defined on O such 
that 

t-l 

Vloga 2 (0)-Vlog5f(0) = ^dipXlogoi^)- logoff)) 

i=l 

P 

+ X)dtj(0)Vlog O ^ 1 _ i (»). 

The sequences of functions (dj), (djj), 1 < J < p, satisfy the same uniform rate of convergence as 
the functions a, Cj+, cx_ and Cjj in (7.7). An application of (7.4) yields to the existence of K > 
and p € (0, 1) such that sup@ ||Vlog<7 2 (0) — V log of (0)|| < Kp l , for almost all trajectories. Point 
3. follows easily and the asymptotic normality is proved. □ 
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8 Conclusion 



In this paper, we investigated the probabilistic properties of the log-GARCH(p, q) model. We found 
sufficient conditions for the existence of moments and log-moments of the strictly stationary solu- 
tions. We analyzed the dependence structure through the leverage effect and the regular variation 
properties, and we compared this structure with that of the EGARCH model. 

As far as the estimation is concerned, it should be emphasized that the log-GARCH model ap- 
pears to be much more tractable than the EGARCH. Indeed, we established the strong consistency 
and the asymptotic normality of the QMLE under mild assumptions. For EGARCH models, such 
properties have only been established for the first-order model and with strong invertibility con- 
straints (see Wintenberger and Cai, 2011). By comparison with standard GARCH, the log-GARCH 
model is not more difficult to handle: on the one hand, the fact that the volatility is not bounded 
below requires an additional log-moment assumption, but on the other hand the parameters are nor 
positively constrained. 

A natural extension of this work, aiming at pursuing the comparison between the two classes 
of models, would rely on statistical tests. By embedding the log-GARCH model in a more general 
framework including the log-GARCH, it should be possible to consider a LM test of the log-GARCH 
null assumption. Another problem of interest would be to check validity of the estimated models. 
We leave these issues for further investigation, viewing the results of this paper as a first step in 
these directions. 
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